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r*** " Abstract 

t— i ; 

To constrain the origin of the soft X-ray excess phenomenon seen in many active galactic nuclei, the 
intensity-correlated spectral analysis, developed by Noda et al. (2011b) for Markarian 509, was applied to 
U \ wide-band (0.5-45 keV) Suzaku data of five representative objects with relatively weak reflection signature. 

They are the typical bare-nucleus type 1 Seyfert Fairall 9, the bright and typical type 1.5 Seyfcrt MCG- 
' 2-58-22, 3C382 which is one of the X-ray brightest broad line radio galaxies, the typical Seyfert-like radio 

loud quasar 4C+74.26, and the X-ray brightest radio quiet quasar MR2251-178. In all of them, soft X- 
ray intensities in energies below 3 keV were tightly correlated with that in 3-10 keV, but with significant 
positive offsets. These offsets, when calculated in finer energy bands, define a stable soft component in 0.5- 
» ■ 3 keV. In each object, this component successfully explained the soft excess above a power-law fit. These 

| components were interpreted in several alternative ways, including a thermal Comptonization component 

i— which is independent of the dominant power-law emission. This interpretation, considered physically most 
reasonable, is discussed from a viewpoint of Multi-Zone Comptonization, which was proposed for the black 
' kJ . hole binary Cygnus X-l (Makishima ct al. 2008). 

■ Key words: galaxies: active - galaxies: individual (Fairall 9, MCG-2-58-22, 3C382, 4C+74.26, 
£f~) , MR2251-178) - galaxies: Seyfert galaxy, Radio galaxy, Radio loud quasar, Radio quiet quasar - X-rays: 

■ galaxies 
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1. Introduction 2011), and the other is a thermal Comptonization compo- 
nent that is separated from that producing the dominant 

In soft X-ray spectra of Active Galactic Nuclei (AGNs), Power-Law (PL) component (Marshall et al. 2003). Since 

a phenomenon called "soft X-ray excess" is often observed, these interpretations often degenerate in spectral analysis 

Characterized by a similar and steep flux upturn towards (Cerruti et al. 2011), there have been no conclusions for 

■ lower energies, this feature is clearly noticed particularly the origin of the soft X-ray excess in AGNs. Recently, 

$_i ' in weakly-absorbed AGNs, such as type I Seyferts and however, Makishima et al. (2008) and Yamada (2011) es- 

Broad Line Radio Galaxies (BLRG). For years, the origin tablished a Multi-Zone Comptonization (MZC) view for 

of this spectral structure has been unidentified, and many Cygnus X-l (hereafter Cyg X-l) through the analysis of 

interpretations have been proposed. Its simplest explana- 0.5-200 keV Suzaku data, and showed that the last inter- 

tion is blackbody radiation from an optically-thick accre- pretation among the three indeed explains the soft excess 

tion disk. However, the observed color temperature of the (Frontera et al. 2001) seen in this leading black hole bi- 

soft excess, typically ~ 0.2 keV, is too high for such disks nary (BHB). This is analogous to the third interpretation 

around black holes (BHs) of > 10 7 M in mass, where Mq of the AGN soft excess as described above, 

is the solar mass. To examine the origin of the soft X-ray excess in AGNs 

Given this, various alternative interpretations were pro- utilizing the wide-band Suzaku capability, Noda et al. 

posed. Among them, mainly three ideas have been (2011b) chose the typical weakly-absorbed type 1 Seyfert 

promising. One is absorption by a partially-covering and galaxy Markarian 509 (Mrk 509 hereafter), and developed 

ionized absorber often incorporating velocity smearing of- a method to study how its 0.5-3 keV intensity in the 

fects (e.g., Middleton & Done 2004; Schurch & Done 2008; Suzaku data is correlated to that in 3-10 keV, and found 

O'Ncil et al. 2007), another is rclativistically blurred re- that a significant positive offset remains in the softer band 

flection from an ionized disk which is in similar condi- intensity when the correlation is extrapolated to lower 

tions as the first case (Zoghbi et al. 2010; Nardini et al. counts. Because the method is based on Count-Count 
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Several AGNs 

Table 1. Information of the five AGNs to be studied. 



Object name 


Type 


Obs. date 


rcdshift 


iV H (Gal)* 


reported Ftf 


previous Suzaku study 


Fairall 9 


Syl 


2007 June 7 
2010 May 19 


0.047 


0.031 


52+ - 20 

-i rc+0.26 
1.00_g 24 


Patrick ct al. (2011) 
Patrick et al. (2011) 


MCG-2-58-22 


Syl.5 


2009 December 2 


0.047 


0.027 


0.69 ±0.05 


Rivers et al. (2011) 


3C382 


BLRG 


2007 April 27 


0.058 


0.074 


0.15±0.05 


Sambruna et al. (2011) 


4C+74.26 


RLQ 


2007 October 25 


0.104 


0.119 


0.3-0.7 


Larsson ct al. (2008) 


MR2251-178 


RQQ 


2009 May 7 


0.064 


0.024 


<0.2 


Gofford ct al. (2011) 



* Equivalent hydrogen column density of the Galactic line-of-sight absorption in 10 22 cm 
' Reflection fraction defined by R = Q,/2ty, when fi is the reflection solid angle. 



Correlation with Positive Offset, hereafter, it is called 
"C3PO" method for simplicity. Thus, utilizing time varia- 
tions with the C3PO method, they extracted a stable soft 
X-ray component, and successfully explained it by a ther- 
mal Comptonization which is independent of the domi- 
nant PL continuum. When this soft Comptonization com- 
ponent is included, the time-averaged 0.5-45 keV spec- 
trum of Mrk 509 was reproduced in terms of a weakly 
absorbed single PL and its reflection. In addition, Noda 
et al. (2011b) discovered that the new soft component 
varied on a timcscalc longer than 3 days, independently 
of the PL variation. This securely excluded a competing 
interpretation of the detected soft component in terms 
of some largely extended thermal emission from the host 
galaxy. 

At the same time as Noda et al. (2011b), Mehdipour et 
al. (2011) analyzed multi- wavelength data of the same ob- 
ject, Mrk 509, which were obtained in a large campaign 
including XMM-Newton, Hubble, and FUSE. They ob- 
tained the same conclusion, that the soft X-ray excess of 
this AGN is created as a thermal Comptonization com- 
ponent which is independent of the principal PL con- 
tinuum. Importantly, the parameters of the thermal 
Comptonization they derived, including an electron tem- 
perature of ~ 0.2 keV, and an optical depth of ~ 16, 
agree with those of Noda et al. (2011b). Thus, through 
the two independent methods (timing analysis and multi- 
wavelength spectral fitting), the soft excess phenomenon 
in Mrk 509 has been confirmed to arise as Comptonization 
by a warm corona, which is presumably different from 
(though possibly related to) the hotter corona produc- 
ing the harder PL continuum. When combined with 
the Suzaku results on Cyg X-l quoted above, this re- 
sult strengthens the general analogy between BHBs and 
AGNs. 

The next step is to examine whether or not the same 
phenomenon as found in Mrk 509 and Cyg X-l is present 
in a larger number of AGNs, including Seyferts and ob- 
jects of other types. For this purpose, we utilize Suzaku 
archival datasets of AGNs, because the wide-band simul- 
taneous coverage with Suzaku, typically available over a 
0.5-45 keV band, is essential in determining the under- 
lying continuum and the reflection component for each 
AGN, and hence to unambiguously identify the soft ex- 
cess signals. In fact, many AGNs have been observed with 



Suzaku since its launch in 2005. However, the overall sam- 
ple is neither complete nor homogeneous. Therefore, we 
choose to conduct the present study via the following two 
steps. The first is to define those AGN types which are 
suitable to our purpose. The second is to select, from the 
Suzaku archive, the best target that represents each of the 
selected classes. 

In the present work, errors refer to ±lcr confidence lim- 
its, except model-fitting parameters in XSPEC for which 
90% error ranges are adopted. 

2. Target Selection 

In order for the C3PO method developed in Noda et 
al. (2011b) to be applicable, the target AGNs should 
satisfy several conditions. First, we should exclude jet 
dominated systems like Blazars. This condition enables 
us to treat pure X-ray emission from a close vicinity of 
the central black hole (BH). Second, absorptions in our 
Galaxy and the host galaxies of AGNs should both be low 
enough to detect soft X-ray signals with significant statis- 
tics. Finally, reflection from accretion disks and other 
surrounding matters should be weak: our C3PO method 
utilizes count rates in the 3-10 keV bands as references, 
so that they should be dominated by the continuum sig- 
nals rather than reflection. Taking these conditions into 
account, we have chosen five AGN types; type 1 Seyfert 
galaxies (Syl), type 1.5 Seyfert galaxies (Syl.5), Broad 
Line Radio Galaxies (BLRGs), Seyfert-like Radio Loud 
Quasars (RLQs), and Radio Quiet Quasars (RQQs). 

As representatives of Sy 1, Sy 1.5, BLRGs, RLQs, 
and RQQs, we chose Fairall 9, MCG-2-58-22, 3C382, 
4C+74.26, and MR2251-178, respectively. Fairall 9, 
3C382, and MR2251-178 are ones of the most typical and 
brightest objects among the currently available Suzaku 
datasets of respective types of AGNs. MCG-2-58-22 was 
selected because of its brightness, moderate absorption, 
and the longest net exposure, among the available Suzaku 
datasets of Syl.5 (Rivers et al. 2011). 4C+74.26 was 
chosen because its X-ray spectrum is similar to those of 
Seyferts, and jet contribution to its 0.5-45 keV signal 
is low (Kataoka et al. 2011). These AGNs were com- 
monly reported to have mild reflection components, e.g., 
by Schmoll et al. (2009), Rivers et al. (2011), Sambruna et 
al. (2011), Larsson et al. (2008) and Gofford et al. (2011), 
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Fig. 1. Background-subtracted and dead-time corrected light curves of the five AGNs measured with XIS FI, in 0.5—3.0 keV 
represented by filled circles and 3—10 keV by open circles, shown with a binning of 10 ks. Open squares represent those of HXD-PIN 
in 15—45 keV, multiplied by a factor of 5 and with the same binning, la errors of the data points of XIS FI are all less than 0.04 
cnt s _1 . Dotted lines show the average 3-10 keV count rates. 



who analyzed the same Suzaku datasets as we utilize. We 
summarize these objects in table 1. The Suzaku archive 
provides two datasets of Fairall 9, obtained in 2007 and 
2010. Among them, we chose the first one for our analy- 
sis, because the object exhibited on this occasion a lower 
reflection fraction and a higher 3-10 keV flux than in the 
other. 

3. Data Reduction 

From the Suzaku (Mitsuda et al. 2007) archive, we 
retrieved the XIS (Koyama et al. 2007) and HXD 
(Takahashi et al. 2007) data of the 5 objects selected 
in section 2. These data, all acquired at the XIS nominal 



position, were prepared via version 2.0 processing in 2005, 
version 2.1 in 2007, and version 2.4 in 2009. 

In the present analysis, the data of XIS 0, 2 and 3 (after 
2006, XIS and 3 only), which use front- illuminated (FI) 
CCD chips, were added and used as XIS FI. The data from 
XIS 1 were not analyzed, since it uses a back-illuminated 
CCD, and has a relatively high and unstable background. 
On-source events of the three or two XIS FI cameras were 
extracted from a circular region of 120" radius centered 
on the source. Background events were accumulated on 
a surrounding annular region of the same camera, with 
the inner and outer radii of 180" and 270", respectively. 
The response matrices and ancillary response files were 
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Fig. 2. The CCPs of Fairall 9. Abscissas gives NXB-subtractcd XIS FI count rate in 3-10 keV, while ordinate gives that in (a) 
0.5-3 keV, (b) 0.5-0.8 keV, (c) 0.8-1.0 keV, (d) 1.0-1.2 keV, (e) 1.2-1.4 keV, (f) 1.4-1.7 keV, (g)l. 7-2.0 keV, (h) 2.0-3.0 keV, and 
(i) 15—45 keV. All data are binned into 10 ks. The error bars represent statistical ±1<t range. The dotted straight line refers to 
equation (1). 



produced by xisrmfgen and xissimarfgen (Ishisaki et 
al. 2007), respectively. 

In a similar way, we prepared events of HXD-PIN. Non 
X-ray Background (NXB) contained in the data was esti- 
mated by analyzing a set of fake events which were cre- 
ated by a standard NXB model (Fukazawa et al. 2009). 
The on-source events and the NXB events were analyzed 
in the same manner, and the latter was subtracted from 
the former. In addition, the contribution from Cosmic 
X-ray Background (CXB; Boldt et al. 1987) was esti- 
mated and also subtracted from the on-source data. This 
was conducted using the HXD-PIN response to diffuse 
sources, assuming the spectral CXB surface brightness 
model determined by HEAO 1 (Grubcr et al. 1999): 
9.0 x 10- 9 (£/3 kcV)-°- 29 cxp(-i;/40 keV) erg cm" 2 s" 1 
str _1 kcV -1 , where E is the photon energy. The esti- 
mated CXB count rate amounts to 5% of the NXB signals, 
in agreement with Kokubun et al. (2007). 

4. Data Analysis and Results 

4-1- Extraction of soft excess component 

While a soft excess component was usually identified as 
excess signals above the PL determined by higher-energy 
spectral data, Noda et al. (2011b) formulated an inde- 
pendent way of doing this. Below, we follow their C3PO 
method. For this purpose, 0.5-3 keV (filled circles), 3-10 



keV (open circles), and 15-45 keV (open squares, from 
HXD-PIN) light curves of the five AGNs are shown in fig- 
ure 1, after subtracting the NXB. All of the XIS FI light 
curves exhibit more than 10% variations, with apparently 
tight correlations between the two XIS bands. To quantify 
this behavior, we produced Count-Count Plots (CCPs), in 
which abscissa is the 3-10 keV band XIS FI count rates, 
and ordinate is those in the 0.5-3 keV. The results, shown 
in panel (a) of figure 2 to figure 7, indeed indicate that 
the 0.5-3 keV count rate, to be denoted y, varies linearly 
with that in 3-10 keV, x. 

For further quantification, we fitted the data in each 
CCP with one straight line, expressed as 

y = ax + b, (1) 

in which the two parameters, a and b, were both left free. 
The fit goodness was evaluated in terms of chi-square 

statistics using errors expressed as a = 



Tl + {aa x ) 2 , 

where a x and a y are 1-sigma statistical errors associated 
with x and y, respectively. 

As shown in table 2 and indicated in panel (a) of figures 
2-6, the data distribution in each CCP has been repro- 
duced successfully with equation (1). Therefore, the main 
variable component is considered to have varied without 
any change in its shape. In addition, the offset b is sig- 
nificantly positive in all cases (2.4-2.8tr), indicating that 
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Fig. 3. The same as figure 3, but of MCG-2-58-22. 
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Fig. 4. The same as figure 3, but of 3C382. 
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Fig. 6. The same as figure 3, but of MR2251-178. Plots without error bars are also shown as insets in panels (b-d). The correlation 
in panel (i) is unconstrained. 
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Table 2. Parameters obtained by fitting 7 CCPs in figure 2(a-i) and figure 3-6(a— h) of the five individual AGNs with equation (1).* 



Object 


Energy range (keV) 


a 


b (counts s 1 ) 


vVd.o.f. 

A. I ^ ■ w • -■- • 


Fairall 9 


0.5-3 


4.36 ±0.10 


0.51 ±0.21 


19.2/18 




0.5-0.8 


0.23 ±0.04 


0.05 ±0.03 


21.1/18 




0.8-1 


0.31 ± 0.04 


0.10 ±0.03 


13.5/18 




1-1.2 


44 ± 05 


06 ± 21 

V7 ■ \J\J — 1 — v ' . — ■ .1. 


16.0/18 




1.2-1.4 


0.37 ±0.06 


0.09 ± 0.05 


22.4/18 




1.4-1.7 


55 ± 07 


02 ± 06 


23 9/18 




1.7-2 


29 ± 04 

V / • — < < — 1 — V / • V 7 I 


05 ± 03 


13.8/18 




2-3 


0.48 ±0.07 


0.06 ±0.05 


24.6/18 




15-45 


0.02 ±0.04 


0.06 ±0.03 


12.3/18 


MCG- 2-58-22 


0.5-3 


2.33 ±0.08 


0.33 ±0.13 


31.5/24 




0.5-0.8 


0.17 ±0.01 


0.01 ±0.02 


40.8/24 




0.8-1 


0.25 ±0.02 


0.08 ±0.03 


39.3/24 




1-1.2 


34 ± 02 


08 ± 03 


44.2/24 




1.2-1.4 


0.37 ±0.02 


0.04 ±0.03 


33.2/24 




1.4-1.7 


46 ± 02 


05 ± 03 


41.2/24 




1.7-2 


29 ± 01 


03 ± 02 


20.8/24 




2-3 


0.49 ±0.02 


0.04 ±0.03 


29.3/24 




15-45 


0.01 ±0.02 


0.15 ±0.03 


36.3/24 


3C382 


0.5-3 


3.78 ±0.26 


0.65 ±0.23 


21.6/23 




0.5-0.8 


0.14 ±0.02 


0.08 ±0.03 


25.4/23 




0.8-1 


0.25 ± 0.02 


08 ± 03 


9.3/23 




1-1.2 


27 ± 03 


17 ± 04 


18.7/23 




1.2-1.4 


38 ± 03 


04 ± 04 

V / . V ' 1 — ! V7 ■ V ' I 


14.1/23 




1.4-1.7 


37 ± 06 

V > • ' J 1 — 1 — V / • V / V7 


18 ± 08 


48 8 /23 

TlV^J wKj I £J *-J 




1.7-2 


27 ±0 03 


06 ± 04 


11 7/11 




2-3 


53 ± 05 


± 0.06 


23 7/23 




15-45 


0.07 ±0.03 


0.02 ±0.04 


25.8/23 


4C+74.26 


0.5-3 


3.22 ±0.30 


0.48 ±0.17 


12.0/14 




0.5-0.8 


0.05 ±0.03 


0.04 ±0.02 


25.6/14 




0.8-1 


0.06 ±0.03 


0.13 ±0.03 


16.6/14 




1-1.2 


0.20 ±0.04 


0.10 ±0.04 


14.8/14 




1.2-1.4 


0.32 ±0.03 


0.04 ±0.03 


8.0/14 




1.4-1.7 


40 ± 07 


08 ± 06 


21.5/14 




1.7-2 


0.19 ±0.05 


0.11 ±0.05 


26.2/14 




2-3 


50 ± 05 


04 ± 05 


11.5/14 




15-45 


0.13 ±0.06 


-0.01 ±0.05 


21.5/14 


MR2251-178 


0.5-3 


1.39 ±0.22 


0.83 ±0.31 


31.4/26 




0.5-0.8 


0.08 ±0.03 


0.07 ±0.05 


19.8/26 




0.8-1 


0.07 ±0.06 


0.17±0.09 


42.9/26 




1-1.2 


0.16 ±0.06 


0.15 ±0.09 


27.7/26 




1.2-1.4 


0.22 ±0.07 


0.11±0.10 


30.2/26 




1.4-1.7 


0.36 ±0.05 


0.05 ±0.08 


12.7/26 




1.7-2 


0.27 ±0.06 


-0.01 ±0.09 


27.4/26 




2-3 


0.26 ±0.10 


0.30 ±0.14 


40.8/26 



* Errors refer to lcr confidence limits. 



some signals should remain in the 0.5-3 keV bands even 
when the intensity of the variable component becomes 
zero. This condition is the same as found with Mrk 509 
by Noda et al. (2011b), and allows us to apply the C3PO 
method to the present five AGNs. 

We divided the 0.5-3 keV band of the individual AGNs 
into seven finer bands, 0.5-0.8 keV, 0.8-1 keV, 1-1.2 keV, 
1.2-1.4 keV, 1.4-1.7 keV, 1.7-2 keV, and 2-3 keV. Then, 
for each AGN, we created seven additional CCPs, where 
abscissa is again the 3-10 keV count rate while ordinate 
is those in the seven finer soft X-ray bands. The obtained 
CCPs are presented in panels (b) to (h) of figure 2 to 6. 



To examine the data distribution in each CCP repre- 
senting a finer energy band, the fitting with one straight 
line expressed by equation (1) was again performed. The 
obtained parameter values are summarized in table 2. 
Thus, the fit with equation (1) was successful in most 
of the 35 CCPs (the seven bands of the five objects). In 
all objects, furthermore, the CCPs generally exhibit pos- 
itive offsets (b > 0), particularly in lower energy bands. 
These offsets, divided by the corresponding energy inter- 
vals, define a spectrum that represents the non- varying 
soft X-ray signals of each AGN. In figure 7, we show these 
spectra in purple, after further normalizing to unabsorbed 
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Table 3. Parameters, with 90%-confidence errors, obtained by fitting the SSE spectra of the five AGNs in figure 7. 



Model 


Parameter 


Fairall 9 


MCG- 2-58-22 


3C382 


4C±74.26 


MR2251-178 


powerlaw 


r 

/V* 


2.89itf 4 

i Q9+0.71 


2 Q2 +i uu 

1 /l +0.60 
A - 4J --0.61 


3 16 +iIb 

°- lu -0.93 
o ftl+0.74 
2 - bi -0.75 


o rw+u.ya 

°- zq: -0.91 
^ 90 4- Q9 


q -1 q + 1.55 
J - ±J -1.80 




x7<io.f. 


2.28/5 


4.09/5 


6.88/5 


7.36/5 


4.50/5 


diskbb 


T disk (keV) 


o.35+s:S 




0.33l^ 13 


0-26to.og 






7V disk 




in r i 7 q 
ly.O zt / .o 


ik Q-i- im 
oo.y zt iu.o 


1 Pi7 1^/1/17 
10 J . 1 zt 44. ( 






2 /J £ 

X /d.o.f. 


3.23/5 


2.74/5 


6.97/5 


7.91/5 


4.88/5 


apec 


T c (keV) 

at! 

/V 4 


i 0+4.8 

3.29t?;£j 


i i+U.5 

lis 

2.78_ 107 


-, 7 +i.o 

1 " ' -0.4 

rr i i 4-2 29 

7.HI2.76 


1 q+0.4 
1 -°-0.6 
7 q-,4-3.24 
'•J'_2.97 


1 -°-0.8 
„ h ,+5 43 

7-1414.25 




x 2 /d.o.f. 


3.49/5 


0.25/5 


7.29/5 


4.48/5 


4.65/5 


kdblur * reflionx 




rn „4-4.y 


50.1+? 7 y 5 


on o + y.o 
60.3_ 7 g 


36. 3+ 7 3 


22.8+g,; 




Aril 

reflionx 


0.62 ±0.23 


0.42 ±0.17 


0.61 ±0.18 


1.97 ±0.47 


5.09 ±1.90 




X 2 /d.o.f. 


2.37/5 


7.38/4 


7.29/4 


4.48/4 


4.65/4 


comptt 


T e (keV) 
r 


0.80tS;3S 
16.2+^ 


0.26i u ;i 7 


0.46'^ 
23.9l 8 7 3 3 


„ „„+u.6U 
u ' '-0.28 

14.9lJ;» 


0.86^;f 3 

1 r 9+179.2 




AT# 

- 1 * comp 


0.21 ±0.08 


0.03 ±0.01 


0.41 ±0.12 


0.83 ±0.24 


0.42 ±0.14 




X 2 /d.o.f. 


2.39/4 


2.59/4 


6.63/4 


7.35/4 


4.57/4 



L" 1-94+Q.45 0.42 ±0.12 3.66+p ^ 6.47+ 2 ; 



* The cutoff power-law normalization at 1 keV, in units of 10 -3 photons keV -1 cm -2 s _1 at 1 keV. 

* The diskbb normalization, in ((_Ri n /km)/(D/10 kpc)) 2 cosf?, where i?i n is the inner disk radius, D the distance to the source, 
and 9 the angle of the disk. 

* The apec normalization, in units of (10~ 17 /47t(Z)a(1 + z)) 2 ) Jn e nndV, where Da is the angular size distance to the source 
(cm), n e and uh are the electron and H densities (cm -3 ). 

§ The ionization parameter, in units of erg cm -2 s" 1 . 

" The reflionx normalization, in units of 10~ J photons keV" 1 cm" 2 s -1 at IkeV. 

* The comptt normalization, in units of photons keV -1 cm -2 s -1 at 1 keV. 
** The 0.5-3 keV luminosity of the extracted SSEs, in units of 10 43 erg s _1 . 



PL models with a photon index T = 2 of which the nor- 
malizations were chosen to approximately reproduce the 
observed time-averaged spectra of the corresponding ob- 
jects. 

To quantify these stable soft emission (hereafter SSE) 
spectra, we fitted them with five representative models; 
power law (powerlaw), multi-color disk (diskbb; Mitsuda 
et al. 1984) emission, thin-thermal plasma radiation 
(apec; Smith et al. 2001), rclativistically smeared and 
ionized reflection (kdblur * reflionx; Laor 1991, Ross 
& Fabian 2005), and thermal Comptonization (comptt; 
Titarchuk 1994). The fits all include the Galactic linc- 
of-sight absorption, which is modeled by wabs (Morrison 
and McCammon 1983) with the TVh values given in ta- 
ble 1. In fitting with the apec model, the abundance 
parameter was fixed at 0.5 solar, while the temperature 
and the normalization were left free. In the kdblur * 
reflionx model, the photon index of the primary contin- 
uum for reflection was fixed at 2, and the inner and outer 
disk radii at 1.24 i? g and 400 R g , respectively. The emis- 
sivity index was also fixed at 4. The other parameters, 
the inner radius, the ionization parameters, and the nor- 
malization, were all left free. In the comptt fit, the seed 
photon temperature was fixed at 0.02 keV, while the other 
parameters were left free. The redshifts in those models 
were fixed at the values given in table 1. The obtained 



fit parameters are summarized in table 3, which also gives 
the 0.5-3 keV luminosity of this component calculated by 
the apec model. Thus, the five SSE spectra can be re- 
produced successfully by any of the employed five models. 
They are characterized by a PL photon index of ~ 3, or 
a disk temperature of ~ 0.3 keV, or an apec tempera- 
ture of ~ 1.2 keV. If employing the relativistically blurred 
and ionized reflection model, the ionization parameter be- 
came ~ 50 erg cm -2 s _1 , with the inner disk radius fixed 
at 1.24 i? g . When employing the thermal Comptonization 
model with a seed photon temperature of 0.02 keV, which 
can be calculated from the mass of a typical central BH 
of ~ 10 8 M Q and an Eddington ratio of ~ 10%, the coro- 
nal temperature becomes 0.3-1 keV and the optical depth 
> 15. Below, wc do not consider the diskbb modeling, 
since its temperature is too high as already pointed out 
previously (section 1). The PL modeling may not be con- 
sidered cither, because the indicated steep PL does not 
accept easy interpretations. In other words, we retain the 
apec, kdblur * reflionx, and comptt modelings. 

In an attempt to extend our analysis to the HXD data, 
we additionally show, in figure 1, 15-45 keV light curves 
of the five AGNs obtained with HXD-PIN. The CCPs of 
these bands against the 3-10 keV XIS counts are shown 
in figure 2— 6(i). The results of fitting these CCPs with 
equation (1) are given in table 2, except for MR2251-178 



No. 



H. Noda et al. 



(b) MCG-2-58-22 




Energy (keV) 

Fig. 7. The High (black) and Low (red) spectra of the present sample AGNs, shown as a ratio to a common unabsorbed T = 2 PL. 
Purple shows the SSE and SHE spectra of each object determined by b of table 2, presented after divided by the same PL as used 
to normalize the High and Low spectra. 



in which the data points have too large errors and the 
variation is too small to derive meaningful constraints on 
the parameters a and b. Thus, the hard-band CCPs of 
Fairall 9 and MCG-2-58-22 suggest statistically significant 
positive offsets (though with large errors), so we divided 
them by the energy range of 20 keV, and show the results 
in figure 7 as ratios to the same r = 2 PL. Hereafter, these 
hard-band components derived with the HXD vs. XIS 
C3PO method are called stable hard emissions (SHEs). 



4-2. High, low, and defference spectra 

In figure 7, time-averaged and background-subtracted 
spectra of our AGNs are also shown, normalized to the 
same (including the normalization) T = 2 PL as used to 
display the SSE and SHE spectra. However, instead of 
summing all photons of each AGN into a single spectrum, 
we accumulated them over two phases, hereafter called 
High and Low phases, when the 3-10 keV source intensity 
is higher (black) and lower (red) than the average (dotted 
line in figure 1), respectively. 
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Fig. 8. Difference spectra of the AGNs in y_F„ forms, obtained by subtracting Low-phase from High-phase spectra, and fitted by 
wabs * cutoff pi. The lower panels show residuals from the best-fit model. 



To examine shapes of the main variable spectral compo- 
nents, presumably PL-like, we subtracted Low from High 
spectra, and show the derived difference spectra of the 
five AGNs in figure 8(a-e). Each difference spectrum 
was fitted successfully by an absorbed cutoff PL model 
(wabs * cutoff pi in XSPEC), in which the photon in- 
dex and normalization, as well as the absorbing column 
TVh , were left free, while the cutoff energy was fixed at 200 
keV. Therefore, in these AGNs, the varying signals can be 
expressed within errors by an absorbed PL (with an as- 
sumed cutoff), in agreement with some previous studies 
(e.g., Schmoll et al. 2009; Rivers et al. 2011; Sambruna 



ct al. 2011; Larsson et al. 2008; and Gofford et al. 
2011). This also explains the good linearity found with our 
CCPs. The obtained parameters are summarized in table 
4. The weakly-absorbed AGNs, Fairall 9, MCG-2-58-22, 
and 3C382, indeed exhibit low absorptions, which mostly 
agree, within their 90% errors, with the Galactic line-of- 
sight column densities reported by Dickey and Lockman 
(1990) and quoted in table 1. Conversely, the strongly- 
absorbed AGN, 4C+74.26, exhibits significantly higher 
absorption than the Galactic value (table 1). 



H. Noda et al. 

Table 4. The parameters obtained in the fits to 0.5—45 keV difference spectra of the five AGNs. 
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Fairall 9 


MCG- 2-58-22 


3C382 


4C+74.26 


MR2251-178 


wabs 




<0.10 


<0.03 


0.08t - 


43 +u - uv 


< 0.50 


powerlaw 


r 




i-8ol ;S? 

0.32+1.1? 


i-63±°:ii 

64+ 017 


i.86t° ;£ 
0.85t°j° 


1 40+ 043 

o ni7+ ' 018 

u - ui ' -0.006 


xVd.o.f. 




58.13/55 


100.15/90 


35.97/55 


35.09/55 


28.43/32 



* Equivalent hydrogen column density in 10 22 cm 2 . 

' The cutoff power-law normalization at 1 keV, in units of 10 -3 photons keV -1 cm -2 s 



4--S. Time-averaged spectrum analysis 

In subsections 4.1 and 4.2, we followed the C3PO 
method of Noda et al. (2011b), and utilized time vari- 
ations to extract the three distinct spectral components; 
the SSE and SHE components (purple in figure 7), and the 
variable (absorbed) PL component (figure 8). Then, can 
we arrive at the same spectral decomposition by analyz- 
ing the time-averaged (i.e., High-phase plus Low-phase) 
wide-band spectra of our sample objects? To see this, 
we fitted these spectra with a model of wabs * zxipcf 
* (cutoffpl + pexrav + zgauss). Here, pexrav rep- 
resents a reflection continuum from neutral matters 
(Magdziarz & Zdziarski 1995), and zgauss an Fel Ka 
line, zxipcf (Reeves et al. 2008) represents the absorp- 
tion by an ionized (or "warm") absorber, with its equiv- 
alent hydrogen column density, the ionization parameter 
£, and covering fraction left as free parameters. This fol- 
lows Weaver et al. (1995), Sambruna et al. (2011), Kaspi 
et al. (2004), and Gofford et al. (2011), who reported 
the presence of warm absorbers in MCG- 2-58-22, 3C382, 
4C+74.26, and MR2251-178, respectively. For Fairall 9, a 
broad iron line component laor (Laor 1991) was added, 
while the zxipcf factor was excluded because this is con- 
sidered to be a bare- nucleus object (e.g., Schmoll et al. 
2008; Patrick et al. 2011). In addition, for MR2251- 
178 which was reported to exhibit absorption by outflows 
(Gofford et al. 2011), we added three gaussians with 
free a and negative normalization, to represent Fei_xvi 
M-shell unresolved transition forest (at a center energy 
of E c = 0.77 keV), as well as those from Fcxxiv L-shcll 
(E c = 1.29 keV) and Fe X xv-xxvi K-shcll (E c = 7.57 keV), 
and fitted with a model of wabs * zxipcf * (cutoffpl 
+ pexrav + zgauss + 3 x zgauss). The cutoff ener- 
gies were fixed at 200 keV, except for MR2251-178 in 
which it was fixed at 100 keV after Gofford et al. (2011). 
The redshifts in these models were all fixed at the val- 
ues given in tabic 1. As a result, we obtained x 2 /d.o.f. 
of 1105.61/658 for Fairall 9, 1099.49/844 for MCG-2-58- 
22, 1070.05/846 for 3C382, 666.18/605 for 4C+74.26, and 
1581.01/871 for MR2251-178: all the five AGNs, except 
for 4C+74.26, rule out the simple PL-based modeling. 
These fit failures in fact arise mainly due to data excess 
above the model in soft X-rays, which does not vanish 
even when the absorbing column density is left free to 
vary. Through this "static" spectral analysis, we have 
thus confirmed the presence of a soft excess component, 
in at least four of our sample objects. 



In order to improve these spectral fits, we tried 
adding a soft excess model, expressed either by a thin- 
thermal plasma model, an ionized and relativistically- 
blurred reflection model, or a thermal Compton model, 
which can represent the SSE (subsection 4.1). Because 
hot plasma emission in the host galaxies is gen- 
erally expected in the soft X-ray band of AGN 
signals, we first employed the thin-thermal plasma 
model, apec, and revised the fitting model to a form 
of wabs * zxipcf * (cutoffpl + pexrav + zgauss + 
apec). As before, Fairall 9 and MR2251-178 ware fitted 
with wabs * (cutoffpl + pexrav + zgauss + laor + 
apec), and wabs * zxipcf * (cutoffpl + pexrav + 
zgauss + 3 x zgauss + apec), respectively. In these 
fits, the temperature, abundance and normalization of 
apec were all left free. Then, as shown in table 5 and 
figure 9, the time-averaged spectra of Fairall 9, MCG-2- 
58-22, 3C382, and 4C+74.26 were reproduced successfully, 
while that of MR2251-178 was not, mainly due to resid- 
uals in the soft X-ray band. Furthermore, even in the 
successful four objects, the obtained abundances (table 5) 
are much lower than the solar value, making the fits phys- 
ically unrealistic. Thus, the thin-thermal emission model, 
which was allowed by the SSE modeling, becomes unsuc- 
cessful when the time-averaged spectra are considered. 

As a next attempt, we replaced the apec model with 
an ionized reflection model, ref lionx, but first without 
relativistic blurring, and tried to reproduce the time- 
averaged spectra with wabs * zxipcf * (cutoffpl + 
pexrav + zgauss + ref lionx), except for Fairall 9 
which was fitted with wabs * (cutoffpl + pexrav + 
zgauss + laor + ref lionx), and MR2251-178 with 
wabs * zxipcf * (cutoffpl + pexrav + zgauss + 
3 x zgauss + ref lionx). The ionization parameter 
and normalization of ref lionx were left free, while the 
photon index of its primary continuum was tied to that 
of cutoffpl. Because the ionized reflection model is 
inevitably accompanied by a strong hard X-ray hump 
even when £ is made very large (e.g., > 10 3 ), the spectral 
soft excess was not fully accounted for by the ref lionx 
model. Specifically, the fits were unsuccessful except for 
4C+74.26, with x 2 /d.o.f. of 1083.54/657, 1049.98/842, 
1206.17/846, 692.57/603, and 1161.51/866 for Fairall 
9, MCG-2-58-22, 3C382, 4C+74.26, and MR2251-178, 
respectively. The fits are thus even worse than those with 
the apec modeling. Given this, the ref lionx component 
in the above fits were convolved with the relativistic 
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Fig. 9. Fits to the time-averaged spectra of the live AGN, with models of wabs * (cutoffpl + pexrav + zgauss + laor + 
apec) for Fairall 9, wabs * zxipcf * (cutoffpl + pexrav + zgauss + 3 X zgauss + apec) for MR2251-178, and wabs * zxipcf 
* (cutoffpl + pexrav + zgauss + apec) for the others, respectively. The SSE and SHE components are also shown in purple. 



smearing kernel, kdblur, with the emissivity index and 
the inner disk radius left free, while the outer radius fixed 
at 400 R g . Then, as summarized in table 6 and shown in 
figure 10, the fits to the time-averaged spectra of all the 
five AGNs were much improved and became acceptable 
for four sources, except MR2251-178. 

In the above successful fits to the 4 sources (figure 10), 
the soft-band (< 2 keV) excess is explained by the kdblur 
* reflionx component, while that in the hard band 
by kdblur * reflionx + pexrav. Furthermore, the PL 
continua of the 4 sources (table 6), against which these 



soft and hard excesses are defined, are found to be con- 
sistent with their High-Low difference spectra (the main 
variable component; tabic 4, figure 8). Therefore, the 
two reflection models, jointly describing the static excess 
above the PL, should also be able to explain their SSE and 
SHE data points which were determined by out dynam- 
ical C3PO analysis. To examine this idea, we compared 
the kdblur * reflionx + pexrav components in figure 
9 with the SSE+SHE data points (in purple), in terms 
of chi-square evaluation but without any parameter read- 
justment. As a result, good agreements were obtained in 
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Table 5. The results of fitting the time-averaged spectra of the five AGNs, wherein the soft excess is represented by an apec model. 



Model Parameter Fairall 9 MCG-2-58-22 3C382 4C+74.26 MR2251-175 



wabs 
zxipcf 


Nn 
Cvr frac. 


0.032(fix) 


0.035(fix) 
7.27l^ 2 
4-35lJ;f 7 
0.581°;^ 


0.074(fix) 
07+ 19 

U.U/ -0.02 
n cc+0.53 
z - oo -0.16 

>0.62 


21 +U U3 
u- zl _o.oi 

1 36+ - 28 

i " TO _0.05 

tS3+ 011 

U -< JO -0.09 

61 + 04 

u - u± -0.05 


0.028(fix) 
1.07 ±0.02 
40+ - 01 
86+ 01 

U.OU_Q Q2 


powerlaw 


r 


1 9Q+ 002 
0.78 ±0.01 


1.85 ±0.01 
1 Q4+0.02 


1.82 ±0.01 
1.11 ±0.01 


1 92+ 08 
j-» z _o.o8 

i in+ - 06 


1.72 ±0.01 

1 I c+0.06 
J- ±o -0.01 


pexrav 


/ref 


i q+o.4 


1.8 ±0.1 


5+ 01 
u.o_ 2 


9 +0 ' 4 

U ' y -0.5 


<0.2 


Fe I Kq 


El 
a (keV) 

iv Fc 

EW (eV) 


6.37 ±0.02 
10- 4 (fix) 
2 30+ ' 21 

Z.OU_ 22 
77+40 
' ' -62 


6.25 ±0.02 
10- 4 (fix) 
1.84 ±0.32 

nn + 54 
OJ -32 


6.41 ±0.03 
0.099 ±0.032 

9 9 r+0.44 
Z - ZO -0.43 


6.40 ±0.04 

10" 4 (fix) 
i 4 2+o.35 

30_? 2 


10- 4 (fix) 

1 ir+0.30 
1 - ld -0.29 

2ii1 


laor 


Rin (Rg) 

TV* 
EW (eV) 


6.37 ±0.07 
16.7lf 8 ° 
1 67+ ' 57 

1 -°'-0.50 

204+^ 










apec 


kT (keV) 
A (Z ) 

»r** 
± * apec 


n o 9 +0.02 
u - oz -0.04 

< 0.005 

^2+°- 15 

U - Jz -0.12 


n 1O+0.02 
U - iS -0.03 

< 0.007 

i - 4i -0.42 


o.25lS:S 

< 0.004 

i ri9 +0 - 45 


0.68±g : « 
< 0.005 

n 70+0. 09 

u - ' 'J-0.12 


0.31 ±0.01 
0.002 ±0.001 
8.06t ; 47 7 


xVd-o.f. 




728.62/651 


846.70/832 


883.80/842 


612.45/601 


983.20/865 



* Equivalent hydrogen column density in 10 22 cm -2 for the Galactic or the intrinsic line-of-sight absorption, 
t The power-law normalization at 1 keV, in units of 10~ 2 photons keV -1 cm -2 s _1 at 1 keV. 

* Center energy in keV in the rest frame. 

^ The Gaussian normalization in units of 10~ 5 photons cm~ 2 s _1 . 
" Center energy in units of keV. 

* The loar Normalizatiation in units of 10 -5 photons cm~ 2 s _1 . 

** The apec normalization, in units of (10 _15 /47t(Da(1 + z)) 2 ) J n c nndV, where Da is the angular size distance to the source 
(cm), n e and riH are the electron and H densities (cm -3 ). 

Fairall 9 and 3C382, with x 2 /d.o.f. of 4.95/8 and 11. 14/8, 
respectively, while poor results in MCG-2-58-22 (19.79/8) 
and 4C+74.26 (23.50/8). Thus, the kdblur * ref lionx 
+ pexrav modeling failed to explain the time-averaged 
spectrum of one source (MR2251-178), and caused incon- 
sistency between the static and dynamic results on two 
more sources (MCG-2-58-22 and 4C+74.26). 

Finally, we replaced the kdblur * reflionx model 
with the thermal Comptonization model, comptt, and 
fitted the time-averaged spectra with a model of wabs * 
zxipcf * (cutoff pi + pexrav + zgauss + comptt), 
except for Fairall 9 which was fitted with wabs * 
(cutoff pi + pexrav + zgauss + laor + comptt), 
and MR2251-178 with wabs * zxipcf * (cutoffpl 
+ pexrav + zgauss + 3 x zgauss + comptt). In 
these fits, the disk temperature was fixed at 0.02 keV 
as before, and the electron temperature of comptt at 
0.5 keV, following Noda et al. (2011b) and Mehdipour 
et al. (2011), while its optical depth and normalization 
were left free. Like in the apec modeling (figure 9), any 
hard-band excess above the PL is now represented solely 
by cold reflection (pexrav). As summarized in table 7 
and shown in figure 11, all the fits are acceptable, and 
the fit goodness is significantly (except 3C382) better 
than those in table 6. Thus, the soft excess component, 



identified statically in each object using its time-averaged 
spectrum, can be reproduced by the Comptonozation 
model. Furthermore in these fits, the comptt component 
has a very similar shape and intensity to the SSE 
spectrum. Actually, when the best-fit comptt model is 
compared (again without any parameter readjustment) 
with the SSE data points, we obtain x 2 /d.o.f. of 4.01/7 
for Fairall 9, 10.02/7 for MCG-2-58-22, 9.52/7 for 3C382, 
9.51/7 for 4CT74.26, and 5.97/7 for MR2251-178, which 
are all acceptable. Therefore, the soft excess structure 
in the time-averaged spectrum and the dynamically 
extracted SSE spectrum can be regarded as two different 
manifestations of the same spectral component, which is 
considered to arise via thermal Comptonization process. 

Similarly, the SHE flux point of four sources (except 
MR2251-178) appear to be consistent with the pexrav 
component. In fact, these SHE data points for Fairall 
9 and MCG-2-58-22 have +2.1 and +5.7 sigma signifi- 
cances above zero, respectively, while they deviate from 
the pexrav model predictions by TO and 2.8 sigmas, re- 
spectively. Similarly, the SHE data points of 3C 382 and 
4C+74.26 (which are both statistically insignificant) de- 
viate from their pexrav predictions by only 0.1 and 0.8 
sigmas, respectively. Thus, the dynamically derived SHE 
can be identified (within the large errors) with the pexrav 
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Fig. 10. Same as figure 9, with the apec model replaced by a kdblur * ref lionx model. The spectrum of Fairall 9 was fitted with 
wabs * (cutoff pi + pexrav + zgauss + kdblur * ref lionx), and that of MR2251-178 was with wabs * zxipcf * (cutoff pi + 
pexrav + zgauss + 3 X zgauss + kdblur * ref lionx) , while those of the others were with wabs * zxipcf * (cutoff pi + pexrav 
+ zgauss + kdblur * ref lionx). 



model required by the time-averaged spectrum. 

5. Discussion 

To study the origin of soft excess phenomena, five 
AGNs with different types were selected from the Suzaku 
archive, and analyzed; Fairall 9 (Syl), MCG-2-58-22 
(Syl.5), 3C382 (BLRG), 4C+74.26 (RLQ), and MR2251- 
178 (RQQ), which have X-ray signals dominated by emis- 



sion from their central engines. Applying the C3PO 
method developed in Noda et al. (2011b) and Churazov, 
Gilfanov, and Revnivtsev (2001), who detected an SSE in 
the bright Syl Mrk 509 and the typical BHB Cyg X-l, re- 
spectively, we succeeded in extracting an SSE in the 0.5-3 
keV band of all the five AGNs. In addition, the applica- 
tion of the same technique to the 15-45 keV HXD-PIN 
data allowed us to detect the SHE from at least two of 
the five AGNs. 
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Table 6. The same as table 5, but the soft excess is represented by ionized and relativistically-blurred reflection (kdblur * ref lionx) 
instead of apec. 



Model Parameter Fairall 9 MCG-2-58-22 3C382 4C+74.26 MR2251-175 



wabs 
zxipcf 


Nn 

N H 

log? 
Cvr frac. 


0.032(fix) 


0.035(fix) 
<0.36 

0.73 + noa 

U.24: 

<0.25 


0.074(fix) 
<0.12 

2.52 ±0.16 
>0.43 


26 +u ' U3 

u - zo -0.02 

24+ 1 11 

u - zy: -0.07 

2.26+°l° 

— U.lo 

> 0.21 


0.028(fix) 
3.15 ±0.02 
2-16+2-S 

— 0.58 

0.65±° ;» 


powerlaw 


r 

iVpL 


1.99 ±0.03 
0.74 ±0.01 


1.84+g-g? 

— U.U1 

1.22 ±0.02 


1.82 ±0.02 
1 04+ 001 


i Q9+0.06 

— U.U4 

n qq+0.08 


1 oq+0.01 

l.OO n no 

1 90+0.04 

1 - ZJ -0.02 


pexrav 




0.8 ±0.3 


< 0.3 


< 0.2 


< 0.5 


< 0.2 


be I Ka 


a (keV) 

N Fc 
EW (eV) 


a qq+0.04 
D -3 8 -0.03 

10- 4 (fix) 
2 45+ - 25 


a on j n no 

6.39 ± 0.03 

10- 4 (fix) 

2.oit2;S 

33±g 


a ^0+0.03 

n95 +0 025 
U.U»O_ 02 g 

9 9 r+0.44 

Z - ZO -0.43 

52^ 


/? on i n n 

6.39 ± 0.04 
10- 4 (fix) 
1 46+ ' 31 

1 - y:O -0.38 

37tg 


c ^n+0-04 
6 - 40 -0.O3 

10- 4 (fix) 
0.91 ±0.25 

otr + 12 
zd -10 


kdblur 


<?* 

-Rin (Rg) 


4.48i ;™ 

9 on+0.49 

^- ou -0.16 


5 15+ 4 19 

i 24+0 35 


2.44t°;j® 


3.96+;-™ 
1 24+ ' 31 


1 24+ ' 27 


ref lionx 


rcnionx 


58.9l3i2 

o i 1+0.34 
• 3 - 11 -0.12 


101.8lgi 

4.97lg;S 


9 oq+o « 


q 4 1+42.9 
»'±- 1 _34.8 

9 44+0.30 


20.2lg 2 2 1 
1 90+0.23 

i - zo -0.09 


x7d.o.f. 




785.06/651 


945.79/832 


886.26/842 


660.51/601 


1145.17/865 



* The emissivity index (scales as R~ q ). 

t The ionization parameter, in units of erg cm -2 s _1 . 

* The reflionx normalization, in units of 10~ 6 photons keV -1 cm -2 s" 1 at 1 keV. 



The highly linear CCPs we obtained rule out possibili- 
ties that the soft X-ray variations of our sample AGNs are 
caused by changes in any absorption, including in particu- 
lar a partially covering warm absorber which is one of the 
three main interpretations explaining the soft excess phe- 
nomena of AGNs (section 1). This is because a partially 
covering absorber, if variable in its covering fraction, col- 
umn density, or ionization degree, would cause complex 
spectral changes, so the CCP would not show such a lin- 
ear distribution as in figure 2-6. Furthermore, in this case 
the soft X-ray count would vanish, without leaving a pos- 
itive offset, when the PL intensity becomes zero. Thus, 
the SSEs must be more independent of the PL continuum. 

The SSE spectra, extracted from the five AGNs via the 
C3PO analysis, were reproduced successfully by five spec- 
tral models; PL, diskbb, apec, kdblur * reflionx and 
comptt. From physical considerations, two of them, the 
PL and diskbb, were ruled out. Combining the static 
analysis with the dynamical one, apec was found to be 
unrealistic. Furthermore, as given in figure 12, the mea- 
sured 0.5-3 keV luminosity of the SSE of our AGNs is too 
high for thermal emission from their host galaxies (which 
would be at most < 10 43 erg s _1 ; e.g., Fukazawa ct al. 
2006), and correlates positively with the AGN luminosity. 
Therefore, the SSE must be tightly connected to the AGN 
phenomenon, rather than to the host galaxies. 

How about the relativistically-smeared and ionized re- 
flection interpretation modeled by kdblur * reflionx? 
When a fine-tuned geometry is actually realized to enable 
the "light bending" effect (e.g., Miniutti ct al. 2007), the 
reflection would not necessarily have to follow variations of 
the primary emission, and may remain stable like the SSE. 



To examine whether or not this interpretation is appropri- 
ate, the broad energy coverage with Suzaku is particularly 
important, because the reflection signals from an ionized 
disk are expected to appear both at the softest and hard- 
est spectral ends. When figure 10 is closely examined, 
we find that the model invoking kdblur * reflionx + 
pexrav tends to undcr-predict the < 3 keV portion of the 
timc-avcraged data, and over-predict the time-averaged 
HXD data. This effect indeed made the fit to MR2251- 
178 unacceptable. Furthermore, at least in MCG-2-58-22 
and 4C±74.26, the kdblur * reflionx + pexrav com- 
ponent determined through our static analysis was not 
able to explain the dynamically determined SSE±SHE sig- 
nals. Thus, the relativistic reflection interpretation can- 
not explain the broad-band Suzaku data of three out of 
the 5 AGNs. In contrast, the thermal Comptonization 
interpretation, using comptt, can consistently explain, in 
all the 5 AGNs, both the dynamically derived SSE and 
the static soft excess. This agreement between the two 
independent methods significantly strengthens the deter- 
mination of the soft excess signals (or the SSE) in our 
sample AGNs. In this case, the SHE can be explained as 
a separate component arising from a cold reflector that 
is most likely located at a considerable distance from the 
central engine. 

Presuming that the SSE is produced as a thermal 
Comptonization component which is separate from (but 
related to) the dominant PL, the corona of each AGN 
is then considered to consist of multiple regions having 
different optical depths and/or temperatures, namely, dif- 
ferent Compton y-parameters. This may be called Multi- 
Zone Comptonization (MZC) condition. Since the various 
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Fig. 11. Same as figure 9, with the apec component replaced by a comptt model. The spectrum of Fairall 9 was fitted with wabs 
* (cutoff pi + pexrav + zgauss + laor + comptt), and that of MR2251-178 was with wabs * zxipcf * (cutoff pi + pexrav + 
zgauss + 3 X zgauss + comptt) , while those of the others were with wabs * zxipcf * (cutoff pi + pexrav + zgauss + comptt). 



types of AGNs which have been selected in the present pa- 
per all have X-ray signals dominated by non-jet continua, 
their central engines (where most of the gravitational en- 
ergy is converted to radiation) are inferred to be in the 
MZC condition. This agrees with the previous results on 
Mrk 509 (Noda et al. 2011b; Mehdipour et al. 2011), and 
the leading BHB Cyg X-l (e.g., Makishima et al. 2008, 
Yamada 2011). The lack of significant short-term (< sev- 
eral hundred ks) variability in the SSE may be explained if 
the SSE-producing Compton corona is largely (> several 



hundred gravitational radii) extended, or more likely, if 
the seed-photon flux (presumably from the disk) is stable 
on these time scales while the dominant PL variations are 
produced in those part of the corona with the largest y- 
paramctcr (Makishima et al. 2008). Overall, we suggest 
that the soft excess emission of some (if not all) AGNs 
is actually a part of its primary continuum, produced in 
grossly the same central engine as the PL component, but 
possibly at somewhat different locations considering the 
clear difference in their variation characteristics. In this 
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Table 7. The same as table 5, but the soft excess is represented by thermal Comptonization component (comptt) instead of apec. 



Model 



Parameter Fairall 9 MCG-2-58-22 3C382 



4C+74.26 MR2251-175 



wabs 
zxipcf 


iV H 
iVn 
log? 
Cvr frac. 


0.032(fix) 


0.035(fix) 
0.36±g;« 
0.36±°;1 4 
0.13+°,^ 

^ ^ — 0.08 


0.074(fix) 
41 + 32 
20+ 16 

°- zu -0.13 
w> 1 -"-—0.48 


0.24+°;- 
u - d ' -0.22 

^ ^ — 0.06 


0.028(fix) 
0.91 ±0.05 
S0+ 04 

o.83 + ?,'!i? 

^ 0.02 


powerlaw 


r 

iVpL 


i s7 +0.04 

0.68±»;° ? 


1 8/L+ 004 

1.33 ±0.01 


1.75 ±0.01 

1.03 ±0.01 


1.84 ±0.06 
0.99 ±0.09 


1.66 ±0.01 

1 04+ 001 


pexrav 


/ref 


0.9+°;! 


1.6 ±0.2 


<0.3 


4+ a4 


< 0.1 


Fe I Ka 


E c 
a (keV) 

iV Fc 
EW (eV) 


fi 07+0. 05 
uo ' -0.04 

10- 4 (fix) 

2.3o±g-.S 

84+" 

— la 


6 26+ 08 

°- zo -0.09 

10^ 4 (fix) 
1.89 ±0.32 


6.43 ±0.03 
0.099l° ;^ 

2.i8tr42 

48tg 


6.41 ±0.04 
10- 4 (fix) 
1 46+ ' 34 
34+1° 


6.40 ±0.04 

10" 4 (fix) 
1 19+0.29 
^^-o.so 


1 3.0 IT 


Rin 

N Fc 
EW (eV) 


u.o/ m u.uu 

16.2+f 9 8 
2 01+ ' 47 

z - ui -0.50 

q-i+51 

yl -65 










comptt 


To (keV) 
kT (keV) 

T 

-^Comp 


-1 7+0.79 
10i ' -0.95 
QO+0.08 


16.2+i; 4 
0.66±° ;i« 


0.02 (fix) 
0.5 (fix) 
1 5 fi+ - 44 

9 97+O.34 
' -0.57 


IS 5+ 2 - 07 

iO -°-1.30 
oc+2.35 


1 7+0.42 
10 - ' -0.54 

42.2+^ 


X7d.cf. 




721.33/652 


825.32/833 


894.07/843 


624.26/602 


929.15/866 



The thermal Comptonization component normalization, in units of photons keV 1 cm s 1 at 1 keV. 
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Fig. 12. The 0.5—3 keV luminosity of the SSE component (table 3) of the five AGNs, plotted against their 3-10 keV luminosity. 
Those of Mrk 509 reported in Noda et al. (2011b) are also calculated and plotted. 



interpretation, the SHE is considered to have a different 
origin from the SSE, and produced via reflection by dis- 
tant cool materials as represented by pexrav. 

The present results have impacts not only on the central 
engine, but also on the determination and interpretation 
of various secondary spectral components, including disk 
reflection, iron lines, and warm absorbers. This is because 



the MZC condition affects the primary continuum shape, 
which was often assumed conventionally as a single PL. 
For example, Noda et al. (2011a) found that the iron Ka- 
line width of the type I Seyfert, MCG-6-30-15, decreases 
considerably, when including a hard component, which 
varied independently of the dominant PL and made the 
primary continuum concave. 
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